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Abstract 

It is well known that in a neutrally-stratified turbulent flow in a deep constant- 
stress layer above a flat surface, the variation of the mean velocity with respect 
to the distance from the surface obeys the logarithmic law (the so-called "law-of- 
the-wall"). More recently, the same logarithmic law has been found also in the 
presence of non flat surfaces. It governs the dynamics of the mean velocity (i.e. all 
the smaller scales are averaged out) and involves renormalized effective parameters. 
Recent numerical simulations analyzed by the authors of the present Letter show 
that a more intrinsic logarithmic shape actually takes place also at smaller scales. 
Such a generalized law-of-the-wall involves effective parameters smoothly depending 
on the position along the underlying topography. Here, we present wind tunnel 
experimental evidence confirming and corroborating this new-found property. New 
results and their physical interpretation are also presented and discussed. 

PACS: 83.10.Ji - 47.27.Nz - 92.60.Fm Boundary layer flows, Near wall turbulence 

In the realm of boundary layer flows over complex topography, much effort has been 
devoted in the last few years to investigate both the detailed form of the surface pressure 
perturbation arising from the interaction between the shear flow and the underlying to- 
pography (see, e.g., Refs. [[[], ||, 0) and its link with the effective parameters describing 
the large (asymptotic) scale dynamics (see, e.g., Refs. @, H). The latter regime is selected 
by observing the flow far enough from the surface and, furthermore, considering solely the 
mean velocity. It is thus clear that with this approach all information on the dynamics 
at smaller scale becomes completely lost. 
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Unlike what happens for the large-scale (asymptotic) dynamics, the description and un- 
derstanding of statistical properties of flows at 'intermediate' scales (a regime which we 
refer to as "pre-asymptotic" , following Ref. ||) seems strongly inadequate. Such regime 
actually attracts much attention in various applicative domains ranging from wind en- 
gineering (e.g., for the safe design and siting of buildings), environmental sciences (e.g., 
for the simulation of air pollution dispersion) and wind energy exploitation (e.g., for the 
selection of areas of enhanced wind speed for the economic siting of wind turbines). 
This almost unexplored regime is the main concern of the present Letter. A first step in 
the understanding of the pre-asymptotic dynamics has been done by the present authors 
in a very recent work ||, where the analysis of simulations of Navier-Stokes flow fields 
[Q over two-dimensional sinusoidal topographies has been performed. More precisely, in 
the case-studies considered, topography takes a sinusoidal modulation of wavelength A 
(along the x— direction, for the sake of simplicity) and amplitude H, its surface having 
an uniform roughness Zq (with z << H). Here, the dominant process governing the 
dynamics is the interaction between the shear flow and the underlying topography, the 
effect of which gives rise to a surface pressure perturbation 0. Such perturbation has a 
depth of the order of H and a downwind phase shift with respect to the topography. The 
latter is the cause of a net force on the flow, acting in the opposite direction of the flow 
itself: thus, an enhanced (with respect to the case of flat terrain) transfer of momentum 
towards the surface takes place. 

Far enough from the surface, the averaged (over the periodicity box of size A) flow will 
'see' an 'effective flat surface' over which the 'basic' logarithmic law (the well known 
"law-of-the-wall" relative to flows over flat terrain |7j]) is restored but now with larger 
(again with respect to the flat case) effective parameters u'f, and zjf ||, on account of the 
enhanced flux of momentum towards the surface originated by the aforesaid shear-flow - 
topography interaction. 

In Ref. 0, we pointed out for the first time - as far as we know - that at least in the 
analyzed WM93 data-set [|J , a generalized law-of-the-wall, is observed: 



where U is the velocity field, x is the horizontal position, z is the height above the terrain, 
and k is the von Karman constant which we will take as 0.4. Notice that the effective 
parameters, wjf (x) and Zq S (x), show a dependence on x at scales of the order of A (i.e. the 
flow 'sees' some details of the topography and not only its total cumulative effects). This 
is precisely the pre-asymptotic regime already defined in Ref. ||. 

In the present Letter, our main goal will be to provide a first experimental assessment 
confirming and corroborating the scenario outlined in Ref. 0. In fact, the main trouble 
of numerical simulations of Navier-Stokes equations is that the impact on the results of 




for 



z > H 
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the closure schemes, through which small scale dynamics is accounted for, cannot be fully 
controlled ||. An experimental confirmation is thus desirable. 

To start our analysis, we briefly describe the experimental set-up relative to the wind 
tunnel experiment performed by Gong et al in Ref. [§]. Details on the description of 
the wind tunnel facility and the basic data acquisition and analysis system are given 
also in Ref. f9j. The experiment was conducted in the AES (Atmospheric Environment 
Service, Toronto, Canada) meteorological wind tunnel, which has a working volume of 
2.44 m x 1.83 m x 18.29 m (wx h x 1). The wave model consisted of sixteen sinusoidal 
waves with wavelength A ~ 610 mm and through-to-crest height H ~ 96.5 mm and was 
placed with its leading edge at distance d ~ 6.1 m downstream from a honeycomb located 
at the downstream end of the contraction region. The topography can be thus considered 
as a fraction of the ideal topography described by: 

h(x,y) = Hsin 2 (^-) (2) 



A 

where y is the perpendicular-to-x-axis direction coordinate. 

Two surface roughnesses were considered, corresponding to the natural foam surface 
(hereafter "smooth case") and to a carpet cover ("rough case"), respectively. For the 
smooth case, velocity profile measurements gave zq ~ 0.03 mm, while for the rough case 
z ~ 0.40 mm. The flow was neutrally stratified and can be considered as a perturbation 
to a stationary horizontally homogeneous infinitely deep unidirectional constant-stress- 
layer flow above a plane surface of uniform roughness, z . Thus, this basic flow should 
have a logarithmic mean velocity profile, U(z) = (ti*/«) \o.{z/zq). The values of were 
~ 0.43 m/s and ~ 0.62 m/s for the smooth and the rough case, respectively. The free- 
stream velocity, Uq, at approximately 1 m above the floor of the tunnel, was set to about 
10 m/s during the measurements both in the smooth and in the rough case. The bound- 
ary layer height, hs, was evaluated to be ~ 600 mm. The rotation of the flow with the 
height, produced in the numerical simulations M by the Coriolis force, is obviously not 
present in the wind tunnel experiment and thus the flow is parallel to the x-axis at all 
elevations. 

Measurements taken over the crests along the hills showed that the flow reached an 
almost periodic state quite rapidly, after the 3rd or 4th wave. Thus, the perturbed velocity 
profiles, U(x, z), were measured at selected downstream locations between the 11th and 
12th wave crests and a very good agreement between the profiles over these two crests 
was confirmed. This topography can be thus considered as a good approximation to a 
two-dimensional topography whose shape is described by Eq. (§) 

In order to compare the numerical simulations analyzed in Ref. |J with the results from 
the wind tunnel experiments here shortly described, we used the same approach as by 



Finardi et al. |T(J and Canepa et al. JTTJ] . Accordingly, noticing that in both cases here 
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considered U ~ 10 m/s, we have kept the speeds (including the friction velocities, w*) 
unchanged, while the wind tunnel lengths (and times) have been multiplied by Xw / Ag ~ 
1000/0.6096 ~ 1640, where Xw and Xq are the wavelengths in the Wood numerical 
simulations and in the Gong experiment, respectively. With this change of scale, hs ~ 
1000 m and A ~ 1000 m in both cases, while the roughness lengths become ~ 0.05 m 
(smooth case) and ~ 0.66 m (rough case), to be compared with the value of ~ 0.16 m of 
the numerical simulations. The hill height, H, becomes ~ 158 m, to be compared with 
the values 20 m, 100 m and 250 m in the Wood numerical experiments. 
The first point to emphasize is that logarithmic laws described by ([I]) are evident also in 
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Figure 1: The local wind speed profiles U(z) from the wind tunnel experiment || are 
plotted (solid lines) as a function of z for four different positions (x in Eqs. (1) and 
(2)) along the hill, corresponding to (a) x — 0, (b) x — A/4, (c) x = A/2 and (d) 
x = 3A/4. The dashed lines represent the unperturbed profile. The dot-dashed lines 
represent the logarithmic law (|l|), with parameters w* ff (x) and z^ s (x) obtained by least- 
square fits performed inside the scaling regions. The values of these effective parameters 
are given in the text. 

the wind tunnel experiments. This can be easily seen in Fig. 1 (the analogous of Fig. 1 in 
Ref. 0), where typical behaviours for the horizontal wind speed profile U(z) (see Eq. (0); 
for the sake of brevity, the dependence on the x-coordinate is omitted in the notation 
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from now on) as a function of z are presented in lin-log coordinates for the rough case 
and for four values of the x-coordinate corresponding to: (a) x = (i.e. h = 0), (b) 
x = A/4 (i.e. h = H/2 upwind), (c) x = A/2 (i.e. h = H) and (d) x = 3A/4 (i.e. h = H/2 
downwind), respectively. Similar behaviours have been found (but not reported here for 
the sake of brevity) also for the smooth case. From this figure, clean logarithmic region of 
the type described by Eq. fll]) are evident and both Wf(x) and Zq C (x) can be measured by 
least-square fits. Specifically, for the four above positions along the hill, we have obtained 
the following values of u^ c (x) and Zq C (x): 1.73 m/s, 0.053 m (for h = 0); 1.38 m/s, 
0.028 m (for h = H/2 upwind); 0.95 m/s, 0.007 m (for h = H) and 1.27 m/s, 0.022 m 
(for h = H/2 downwind), respectively. Such values can be compared with those in the 
absence of any hill: u± ~ 0.62 m/s and zq ~ 0.66 m. 
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Figure 2: The measured (circles) effective parameters u'f(x)/u^ — 1 (on the left) and 
hi(zQ a (x)/zo) (on the right) as a function of the ratio x/X along the axis of the hill in the 
smooth case (above) and in the rough case (below). The continuous lines represent the 
sinusoidal law best-fitting the experimental data. 

The results of the least-square fits are summarized in Fig. 2 where both profiles 
u< fi. x )l u * ~ 1 ( on ^e left) and \n.(zf \x) j 'z ) (on the right) are shown as a function of 
x/X for both the smooth (above) and the rough (below) case (different scales in the or- 



5 



dinates have been adopted). Notice that both uf(x)/u+ — 1 and lia^z^^x)/ z ) have been 
fitted with the analytical expression (0), relative to the topographic profile, but with a 
shift of A/2 (i.e. x i— > x + A/2 in (||)). More precisely, we suggest the expression: 

y(x)= (y) - A y [h(x)/H - 1/2] (3) 

where y(x) stays for either uf(x)/u+ — 1 or In (z^ \x) / z ) , (y) is the average value of y(x) 
in the interval (0, A) and h(x) = h(x,y) is the topographic shape given in Eq. (Q). It 
should be also stressed that we have considered ln(^ ff (x)/zo) instead of the simpler ratio 
Zq S (x)/zq, because the former parameter is more similar to the topography shape than 
the second one. 

It is now interesting to put together the new results here obtained with those of Ref . . 
This allows to investigate at which degree of accuracy one can express the behaviours 
of quantities like the average values, (y), and amplitudes, A y , of these sinusoidal shapes 
solely in terms of simple geometrical parameters. Simple considerations suggest to look 
at the ratio H/ A: this is indeed a rough measure of the hill slope. 
The values of (y) are reported in Fig. 3 for both the smooth and the rough case. The 
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Figure 3: The mean values (u^(x)/Ui)—l and (\n(z'jf (x) / z )) in Eq. (3) versus H/ A. Stars 
are relative to the Wood and Mason numerical simulations (u* ~ 0.44 m/s, zq ~ 0.16 m); 
diamonds are relative to the smooth case («* ~ 0.43 m/s, zq ~ 0.05 m) and squares are 
relative to the rough case (u* ~ 0.62 m/s, zq ~ 0.66 m). Dashed lines represent the linear 
curve best-fitting the Wood and Mason numerical data. 

dashed lines (a linear fit in H/X) are obtained considering only the results of numerical 
simulations, while the values from the wind tunnel experiments are reported with their 
error bars. These monotonic behaviours are expected on account of the increasing of 
the (total) transfer of momentum towards the surface arising for increasing slopes. The 
amplitudes A y are shown in Fig. 4. The dashed lines are a parabolic fit in H/ A and, as for 
(y), they have been obtained by only considering the results of the numerical simulations. 
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The values from the wind tunnel experiments are again presented in the same figure. 
Curves relative to the amplitudes reach a maximum for H/X ~ 0.20, after that start to 
decrease. We can argue that two different mechanisms exist and act in competition. The 
physical key role is played by curvature effects [12|, already invoked in Ref. |J to explain 
the presence of minima (maxima) located above the hill top (valley) for both u'flx) and 
hiz^ i {x). To be more specific, let us consider the two opposite limits H/X <C 1 and 
H/X 3> 1, from which we can easily isolate the two competing mechanisms. Concerning 
the former limit, we have gentle slopes and it is well known that in this case the flow 
closely follows the surface contour. Streamlines are (weakly) curved and, as pointed out 
in Ref. []12 |, energy is transferred towards the large scale components above the hill tops, 



while it blows towards the smaller scales above the valleys. The quantity u'f 2 (x), that 
is a measure of the energy of turbulence, is thus smaller on the hill top than above the 
valley. Let us increase (just a little bit) H/X. The flow again closely follows the surface 
contour but streamlines are now more curved. As pointed out in Ref. [12| , energy transfer 
thus increases and, as an immediate consequence, the same happens for the difference 
between the maximum and the minimum of u"f(x). But this means an augmentation of 
its modulation amplitude. 

In the second limit H/X ^> 1, a further important effect arises due to trapping regions 
placed on the downstream hill slopes. It is in fact well known (see, e.g., [|K|) that, for 
surface slopes large enough, the flow is not able to follow the contour surface and sepa- 
rates. In this case, in between two hill crests, the flow is essentially trapped and, roughly 
speaking, streamlines are expunged in the wake region. The dynamical consequence is 
that the flow streamlines are weakly modulated, and this also happens for the shape of 
u"f(x). If we now decrease the ratio H/X, trapping effects reduce and this means that the 
wake can penetrate more deeply in the valley, with the consequent increasing of curvature 
effects and thus of the u"f(x) amplitude. 

From the inspection of these two limits, it is thus clear that a maximum in the amplitude 
should be attained for a certain finite value of H/X, i.e. when the two competing mecha- 
nisms are balanced. 

Being the maximum (minimum) of In Zq S (x) directly related to the presence of the maxi- 
mum (minimum) of u'flx) (see Ref. || for the discussion of this point) the argumentations 
above presented hold also for \iizq B (x). 

Comparing the values of both amplitudes and mean values of the effective parameters 
extrapolated from the numerical simulations (i.e. from the linear fits in Figs. 3 and 4) and 
those from the wind tunnel experiments, we notice that, for the smooth case, experiments 
are always compatible (within the error bars) with the numerical simulations. This is 
not always the case for the rough case. We remark that u* and zq (relative to the flat 
terrain) are closer to the WM93 case studies in the smooth case than in the rough case. 
This suggest that the expression of both amplitudes and mean values solely in terms of 
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geometrical quantities like the ratio H/ A is a reasonable approximation for small variations 
of the 'bare' parameters and z . When the range of variability of the latter two 
parameters increases, an explicit dependence on them has to be taken into account. 
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Figure 4: The values of the amplitude A y in Eq. (3) relative to y — u"f(x)/u± — 1 and 
y = ln(z5 fi (x) I ' Zq) versus H/X. Stars are relative to the Wood and Mason numerical 
simulations (u* ~ 0.44 m/s, Zq ~ 0.16 m), diamonds are relative to the smooth case 
{u± ~ 0.43 m/s, zo ~ 0.05 m) and squares are relative to the rough case (u* ~ 0.62 m/s, 
z ~ 0.66 m). 

For a better evaluation and understanding of the dependence of (y) and A y on H/X, 
zq and u*, the analysis of more numerical and wind tunnel experiments, and possibly in 
nature, is necessary. Nevertheless, the wind tunnel data here considered give a strong 
confirmation of the existence of a pre-asymptotic regime characterized by a generalized 
law-of-the-wall given by Eq. ([]]) and pointed out for the first time in Ref. 0. Thus, this 
phenomenon appears as a real physical property and not a spurious feature produced by 
some of the approximations (e.g. parameterizations of small-scale, unresolved dynamics) 
used to solve the Navier-Stokes equations. 
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